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■ In this work we derive evolution equations for the nonlinear behavior of a coasting beam under 

' the influence of a resonator impedance. Using a renormalization group approach we find a set of 

coupled nonlinear equations for the beam density and resonator voltage. Under certain conditions, 
these may be analytically solved yielding solitary wave behavior, even in the presence of significant 
dissipation in the resonator. We find long-lived perturbations, i.e. droplets, which separate from 
the beam and decelerate toward a quasi-steady state, in good agreement with simulation results. 
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I. INTRODUCTION. 



o 
o 

CJ 

Observations of long-lived wave phenomena have been made in stored high-energy beams for many years. For 
the most part, these have been ignored or avoided as pathological conditions that degraded the performance of the 
machine. However, in recent experiments, as well as in simulations, observations have been made which suggest the 
£L|) occurrence of solitary waves in high-energy stored beams under certain conditions. Both from the point of view of 
scientific curiosity as well as the importance of understanding the formation of halo in such beams, it is worthwhile 
I ' to study the physics of these nonlinear waves. 
, Of particular interest is the saturated state associated with high-intensity beams under the influence of wakefields, or 
CNj 1 in the frequency domain, machine impedance. In stored beams, especially hadron beams where damping mechanisms 
are relatively weak, a tenuous equilibrium may develop between beam heating due to wake-driven fluctuations and 
damping from a variety of sources. This state may well be highly nonlinear and may depend on the interaction of 
nonlinear waves in order to determine the final equilibrium state. It is our interest in this work to elucidate the 
conditions under which nonlinear waves may occur on a high-energy stored beam. This will then lay the groundwork 
' for a future study of the evolution of the beam under the influence of these nonlinear interactions. 

We note that much work has been carried out already on solitary waves, j5|, ||, [jrj, and references contained 
' therein, including those occurring on a beam under the influence of internal space charge forces Q|, Our situation 
is new in that we consider the specific form of a wakefield associated with a high-energy beam, namely when space 
charge forces are negligible. This leads to a specific form of a solitary wave in a dissipative system, one which has 
received limited attention in the literature thus far ^ , jl(J , [|ll| . We have made both experimental observations and 
carried out simulations which show the long-lived behavior of the nonlinear waves even in this dissipative case. It is 
^ ' our aim to shed light on this case. 

In this work we adopt an approach which is commonly employed in fluid dynamics to arrive at a set of model 
equations for solitary waves on a coasting beam under the influence of wakefields. It is based on the renormalization 
group (RG) analytical approach, which is akin to an envelope analysis of the wave phenomena. The method in the 
form we will use it was introduced by Goldenfeld || and expanded upon by Kunihiro Q . 

In Section II we derive the amplitude equations for a resonator impedance following the standard renormalization 
group approach. This results in a nonlinear set of equations for the wave amplitude and beam density. In Section 
III we proceed to find analytic solutions for this set which does indeed admit solitary waves. In Section IV we give 
the conclusions of this study and outline the procedure for applying these results to the study of the steady-state 
fluctuations on a stored beam. 
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II. DERIVATION OF THE AMPLITUDE EQUATIONS. 



Our starting point is the system of equations 

df df Rf „ 

^+ 2 7^+<^ = ^ (1) 
J(0;T) = / dvvf(9,v;T) 

for the longitudinal distribution function / (0, v; T) of an unbunched beam and the voltage variation per turn V (9; T). 
To write down the equations (|l|) the following dimensionless variables 

rp i 9 k a e uj r w 

T = uj s t ; v = — = H ; w= — ; 7 = (2) 

_ e 2 Rk ~/ 

have been used, where u s is the angular revolution frequency of the synchronous particle, e is the energy error, ujr is 
the resonator frequency, Q is the quality factor of the resonator and R is the resonator shunt impedance. Furthermore 

is the proportionality constant between the frequency deviation and energy deviation of a non synchronous particle 
with respect to the synchronous one, while r\ = olm — % 2 (cum - momentum compaction factor) is the phase slip 
coefficient. The voltage variation per turn V and the beam current / entering eqs. (|l]) have been rescaled as well 
from their actual values V a and I a according to the relations 

V a = 2euj sl RV ; I a = euj. (4) 

Let us now pass to the hydrodynamic description of the longitudinal beam motion 

dp 9 . . 
df + d9 {pu)= °> 

du du o-^dp 
d 2 V n dV 2rr d . , 

where 



p (9; T)= I dvf (9, v; T) ; p (9; T) u (9; T) = J dvvf (9, v; T) (5) 

a v = (6) 

UJ S 

and cr e is the r.m.s. of the energy error that is proportional to the longitudinal beam temperature. Rescaling further 
the variables p and V according to 



2 



Pa = PoP 



V a = 2euj s p jRV 



X 



e 2 R-fk p 



(7) 



and taking onto account that the dependence of all hydrodynamic variables on 9 is slow (~ eff) compared to the 
dependence on time we write the gas-dynamic equations as 

dp d . . 

- + e- ( P u) = 0, 



du du a 2 dp 

df + £u de =xv - £ JW 



(8) 



d 2 V n dV 2Tr d . . 
— +2j—+u J V = —( P u). 

Here e is a formal perturbation parameter, which is set to unity at the end of the calculations and should not be 
confused with the energy error variable. We will derive slow motion equations from the system (||) by means of the 
renormalization group (RG) approach |], f§. To do so we perform a naive perturbation expansion 

oo 

p = 1 + ]T e m Pm ; 

m— 1 

around the stationary solution 



e u 



m i 



V=J2 e m V m 



The first order equations are 
dpi 



dT 



~df~ 



AVi 



V<® = 0. 



d 2 V x n dV x - dm 

+ 2t — - + iJ-V-x = — - 

8T 2 1 dT 1 dT 



with obvious solution 



V 1 {9;T)=E{9;T )e 



E* (0;T o )e 



r at 



(9) 



(10) 



(11) 



ui (9;T) = u (9;T ) + \ 



E_(6jTo)_ e , Wl AT _ ^* ( 6, ; T °) C -^TAT 



(12) 



Pi (6;T) = R (6;T ). 

In expressions dll]-[l3]) the following notations have been introduced 



U>1 = UJ q + 17 



2 2 \ 

uj q =uj — A, 



(13) 



(14) 



AT = T - T n 



(15) 



where the amplitudes Tv (0; T G ), it (9; T D ), i? D (0; T Q ) are yet unknown functions of 9 and the initial instant of time T a . 
Proceeding further we write down the second order equations 

d P2 d 

W + d9 {Pl+Ul)= °> 



du 2 d Ul _ 

~dT ~&T ~ 2 ~ v d9 ' 



2 d Pl 



d 2 V 2 n dV 2 2t , d , 

~df 2 + 1 ~&T + U 2 = dT ^ 2 + PlUl + P2) ' 
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Solving the equation for the voltage 

d 2 V 2 „ dV 2 



that can be obtained by combining the second order equations, and subsequently the other two equations for u 2 and 
p 2 we find 



V 2 (9: T) = - Jj [2u' a + (*l + 1) R>] + ^ ( Ro E + ?tf) - r. 



(16) 



u 2 (6: T) = -L' + alK + A [ 2 < + (a 2 + l) #J ) AT + AtfV- at + 



+ * ( RoE + *^ e — _ ^ f R „ £ + 2* ] e — + c . c . 



2iu q u){ 



(17) 



p 2 (0; T) = — (R' a + u' ) AT + ±E'e^ AT + c.c. 



(18) 



where the prime implies differentiation with respect to 9. In a similar way we obtain the third order equations 

-^+Q- e {p2+RoU l +U 2 ) = 0, 



du 3 dm du 2 2 / dp 2 dR " 



d 2 V 3 . „ 0V 3 



+ 2 7-7^l + u 2 V 3 = {u 3 + R Q u 2 + p 2 ui + p 3 ) . 



d 



dT 2 ' dT 
Solving the equation for the voltage 



dT 



d 2 V 3 , „ dv, 



gji + 2 7 ^f + u 2 V 3 = -2u> 2 - 2u lU [ - (a 2 + l) p' 2 + XR V 2 - 2 (R oUl )' + Xp 2 V 1 

that can be obtained by combining the third order equations, and subsequently the other two equations for u 3 and 
p 3 we obtain 

V 3 (9; T) = - A f 2 u u' + A [2R oU ' + (<j 2 v + l) R R' ] + 2 (R u )'\ 

^ O I ) 



" + alR" + — (2< + (<j 2 + 1) R" Q ) 



(a 2 + l) « + <) 



AT 



>> + alK + - 2 {2u:+(al + l) K) 



+ (a 2 v + 1) « + <) 



+ A 2 



T'T* ST* 



, i* 2 



e -2 7 AT + 



AAT 

2iuo a 



I 9 „\ -E"' A 



i (RoE)' 



2E" 



+ -[K + i? )£]'[e^ AT 
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<^AT , 



LOiLOg 



(RoE)' + 



2iE" 



(jji 



2iuj a 



(u' + R' )E^e^ AT + c.c 



(19) 



U3 (0; T) = (alR R' u u' ) AT + ^ °ff e~^ AT - 
- A | 2u X + ^ [2«' + (a 2 + 1) K] + 2 (iW) AT- 



2 7 A 



^ + aX + 4( 2 <+K + l)<) 



+ (a 2 + l) « + <)} AT 



A | 2A^|^ +A2 ^ + ^^ e _ 27AT+ 



+ oscillating terms and terms proportional to (AT) Z 



(20) 



p 3 (6;T) = -(i? oUo )'T+ 



oscillating terms and terms proportional to (AT)* 



(21) 



Collecting most singular terms that would contribute to the amplitude equations when applying the RG procedure, 
and setting e = 1 we write down the following expressions for Vrg, urg and Prg 

V rg (0;T,T o ) = Ee l ^ AT + 



R E+ — [(1 + u + R ) E] - (al + 3 N 



2iuj a 



U>1 



Uj\ 2uJ q 



(u + R )' E > e iLJlAT 



+ 



A 2 AT f wi + 2lo q 



2iu 2 \ 2ujf 



i (RoE)' - 



2E" 



4uj„ 



(22) 



urg (0;T) 



no - {< + a 2 v R' + A [ 2u ' + (a 2 + l) R' a ] J AT- 



+ K 2 i? < - uX) at + ^^- 27AT - 

~ \2u u' + ^ [2< + (a 2 + 1) + 2 (i? u )'} 



27AAT f 



< + <tX + 4(2< + K 2 + i)K0 



+ (a 2 + l) (u'l + K) 
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3 -2 7 AT 



(23) 



Prg 1 



; T) = R a — [R' Q 



{R u )'] AT. 



(24) 



The amplitudes E, u and R can be renormalized so as to remove the secular terms in the above expressions ( f22f - 
p4| ) and thus obtain the corresponding RG equations. Not entering into details let us briefly state the basic features 
of the RG approach M. The perturbative solution ( p2"f]24| ) can be regarded as a parameterization of a 3D family of 
curves {3?t„} = {Ro (T ) , u (T ) , E (T )) with T D being a free parameter. It can be shown that the RG equations 
are precisely the envelope equations for the one -parameter family {5Rt } : 



dR du Q dE 

#zy Wo 1 df 



o. 



T a =T 



It is straightforward now to write down the RG equations in our case as follows: 

OR d 



(25) 



(26) 



di 



d 



dR n 



~dT + d9^ Uo + ° vRo > + U °~d6 ~ avRo 'de 



A 2 d \E\ 
w 2 36 



-2-yT 



2A 



R' 

u' a + u u' + {R u )' + (cr 2 + l) -± 



\ 2 R. 



°-[2u> +(ol + l) R' ] 



Uu'l + 2*tK + ~ 2 K' + K + 1) K] + {ol + 1) « + <) 



2X 2 d\E\ 2 ~(E'E* 
A 



00 



EE* 



-2-yT 



(27) 



2iu>. 



d d 



2% 



[(1 



Ro)E]' 



-E"- 



2u>, 



{l + Uo + Ro)' E- 



X_ 



Ul + 2u) q 

2uj\(jj q 



i (R E)' 



2E" 



R 



r e 



2iE' 



(28) 



In deriving eq. (|2q ) we have assumed that the voltage envelope function E depends on its arguments as E (6 — T a ; T a ). 
Neglecting higher order terms we finally obtain the desired equations governing the evolution of the amplitudes 



dp 
df 



lfe {pu) 



= 0, 



(29) 



du 
df 



,du 
l d6 



a^dp 
' p d9 



\2 d 



E 
~dT 



(30) 



G 



2ioj q ( d . d . \ ~ ,„ ^ ~ o-, 2 + 3 5 2 £; 



«^ + ^+7 S=(p-1)B- 



A \dT 86 'J VP ; cof d6 2 

2i d / ~\ % -dp 

where 

p=l + R Q ; u = l + u a ; E = Ee~ lT . 
Eliminating u from equations (|2|) and ( |30|) we get 

~ 2 

9 2 p 2 d 2 p \ 2 d 2 E 



— a. 



dT 2 v d6 2 lu 2 d8 2 ' 
2iu g ( d d \~ a 2 + 3d 2 E 2i dE „ ~ 

III. SOLUTION OF THE AMPLITUDE EQUATIONS. 

Let us perform a scaling of variables in the amplitude equations (|3^), (|34| ) according to 

AT w o |A|£ 

t = - — ; 6 = — ; V = • 

2w 9 yV 2 + 3 cr„w 



The amplitude equations take now the form 



d 2 p 2 d^i_ 2 d 2 \^ 



where 



Q T 2 uq q2 u g Q2 



. d abuj d \ u? d 2 ip , . w dip 

i ywT + —d^)^ +iib ^ = -^dQ 2+ia ^d@ +{p - i) ^ 



2 2uj„ 2cr v uj a Lu 
a = — , ; 6 = ■ " — ' 



^f+3 A |A|y^T3' 

From equation (^6|) one finds approximately 

p = i-m 2 . 

Equation ( |39|) when substituted into @) yields 



d abuj d \ a; 2 9 V . w o 



Noting that 



T 

cji = uj e l ^** E ; £Jar g = arctan — 

and introducing the new variable 
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a; = B + ar — r = 



- T 



2 ' 

we rewrite the nonlinear Schrodinger equation (40) in the form [|Io| 

2i7~\ 9 2, 07 dtp 



XT 

LU Q UJ q 



I— h = 

ar 



cj J 9x 2 w dx ^ ^ 



Next we examine the linear stability of the solution 

4> (x-T)=A e l( * kx - nT \ 

where 

n = k 2 - A 2 - — (2k 2 -ak + uj b) . 
In the case the energy of the beam is above transition energy (k a < 0) the solution 



i-Ji+^<^<iWi+^ 



a a 

To proceed further let us represent the field envelope function tp as 

i)(x;r) = A(x;t) e^'^ 
and write the equations for the amplitude A and the phase ip 

A T + jbA = -A<p xx ~ 2A x ip x + — (A xx - Alp 2 x ) + —A<p x , 

Wo Wo 



(41) 



(42) 



(43) 



is exponentially decaying for 
(44) 

(45) 

(46) 



Aip T = A xx - Atp 2 x + A 3 + — (Aip xx + 2A x tp x ) 



a-/ 



A x . 



When 7 = the above system admits a simple one-soliton solution of the form 

ip (x] t) — kx — fir + a, 



(47) 



(48) 



A (2; r) 



V2K 



cosh [K (x - 2kr + 0)] 



K = k — Q. > 0. 



Define now the quantities 



N(t) = / dx\i/j(x;T)\' 



P(T) = 2 ^V-dx"^ dx- 



(49) 



(50) 



These are the first two (particle density and momentum respectively) from the infinite hierarchy of integrals of motion 
for the undamped (7 = 0) nonlinear Schrodinger equation ||]. When damping is present (7 7^ 0) they are no longer 
integrals of motion and their dynamics is governed by the equations 



^ + 2 7 6iV=-^ / dx 
clt uj. 



dtp 



dx 



2«7 



(51) 



dP 
d7 



+ 2jbP 



2ij f (d 2 tpdip* d 2 tp* dip\ 2a 7 



dx 



\ dx 2 dx dx 2 dx 



+ 



dx 



dtp 



dx 



(52) 



Instead of solving equations ([l6]) and (^7|) for the amplitude A and the phase p we approximate the solution of the 
nonlinear Schrodinger equation (E2J) with a one-soliton travelling wave 
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where 

H{t)=2 J 'dra(r) ; Q (r) = J dr [a 2 (r) - if (r)] . 

Substituting the sample solution (Q) into the balance equations (pl|), ( p2[ ) and noting that 

N (r) = 477 (r) ; P (r) = 4tj (r) cr (r) 
we obtain the following system of equations 



dr) 



47 / 77' 



uj V 3 



- — h 2767; = — + rj<7 H 770-, 



207 



lJ 



dr 



27677(7 = (rfa + 7/cr 3 ) H ^- + 7/cr 



or 



d?7 
dr 7 



27677 = — + rja I H 770-, 



^ = _^2L 2 cr+ 2 
dr 3u 3lu 

In order to solve equations (^5h and (|56|) we introduce the new variables 



Ht)=V 2 (r) 



k(t) =ct(t) 



so that the system fl55|), ( p6[) is cast into the form 

3uj oj 

where 



47^-^ 2 -^ 2 



d« 87 
dr 3oj 



h = ^-b>0. 
8uj 

A particular solution of the system of equations ( ^8| ) can be obtained for n — 0. Thus 

?7 2 (0)e 4 T blT 



a 

u — — 
4 



77 (t) = 3W &1 



3w Q 6i + 2t7 2 (0)(e 4 ^ T - 1)' 



Solving equation (|29|) for u, provided p is given by fl39j ) and (53) one finds 

A^/cr 2 + 3 cosh 2 z 



u (x; t) 



2iv iv q (cosh z — 2?7 2 ) 



\~, , fi 2 ?7 2 \ , I6777 3 ,3 i] 2 - co 

C + 4777 61 — tanh z H — tanlr z + a- ; 

V w o / 3w cosh" 



cosh 2 z 



where 



z (x; r) = 77 (r) [x - /i (r) + f3] , 



C = a [1 - V 2 (0)] + 



2aj u; q 



[1 - 2 V 2 (0)] [1 + u (0)] . 



(63) 



The solutions for the mean velocity of the soliton and the corresponding voltage amplitude are shown in Figs. 1 
and 2 respectively. We note that the solitary wave corresponds to a self-contained droplet of charge which separates 
(decelerates) from the core of the beam and approaches a fixed separation at sufficiently long times. The reason for 
this behavior is the fact that the driving force due to the wake decays rapidly as the soliton detunes from the resonator 
frequency. At sufficient detuning, the wake no longer contains enough dissipation to cause further deceleration. The 
resonator voltage decreases in a corresponding fashion. It is interesting to note that the charge contained in the 
soliton remains self-organized over very long times despite the presence of dissipation. This situation is rather unique 
and is due to the peculiar character of the wake force from the resonator. 




T (time) 



FIG. 1. Mean velocity of the solitary wave due to a resonator impedance. Solitons decelerate at first due to the dissipative 
part of the wakefield. However, over long times, they approach a steady state where the wakefields have sufficiently decayed 
due to the finite resonator bandwidth. 
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T (time) 



FIG. 2. Voltage amplitude on the resonator. The voltage first grows due to the longitudinal impedance, followed by 
oscillations which result from the interference of energy between the solitary waves and the core of the beam. The envelope of 
the amplitude eventually decays as detuning occurs. 

In Figs. 3 and 4 we show the corresponding mean velocity and voltage from a coasting beam simulation previously 
reported. The behavior is manifestly similar to that predicted by Eq. (35), (53) and (61) , though no attempt has 
been made to check the precise scaling of the physical quantities. 




-BO »~ 
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lHOOO lSOiMi J0OO0 25000 30000 35000 40000 4500* 5'JGOO 



FIG. 3. Mean velocity of the solitary waves from the simulation showing deceleration toward a fixed maximum energy 
separation. There is good qualitative agreement with the analytical result. 
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FIG. 4. Voltage amplitude on the resonator from the particle simulation. There is good qualitative agreement between the 
analytical results and the voltage envelope shown. 



IV. CONCLUSIONS 



In this work we have derived a set of equations for solitary waves on a coasting beam using a renormalization group 
approach. This procedure has led to a specific set of evolution equations in the practical case of a cavity resonator 
of finite Q. The resulting set of equations can be solved analytically under certain assumptions, and this leads to an 
explicit form for the soliton and its behavior over time. We find, in contrast to other solitary waves in the presence 
of dissipation, that solitons can persist over long times and do so by decelerating from the core of the beam. This 
deceleration leads to detuning and the decay of the driving voltage. The result is that a nearly steady state is reached, 
albeit with a gradually decreasing soliton strength, but fixed maximum energy separation. 

Good qualitative agreement between the analytic results and the simulations have been observed. We note that 
such a process may well indicate a method by which well-defined droplets can occur in the halo of intense stored 
beams. Further study of this problem, and the application of the RG approach to bunched-beam evolution will be 
considered in future work. 
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